Journal of Research of the National Bureau of Standards Vol. 59, No. 3, September 1957 Research Paper 2784 

Estimation of the Frequencies of Thin Elastic Plates 

With Free Edges 

Tosio Kato, 1 Hiroshi Fujita, 2 Yoshimoto Nakata, 2 and Morris Newman 

A variational method is proposed for calculating the frequencies of thin elastic plates 
with free edges, with rigorous error estimates. As a numerical example, the fundamental 
frequency of a square plate with the Poisson ratio 0.225 is calculated with a satisfactory 
result, the possible relative error being less than 1/2000. Generalization to more complicated 
boundary conditions is straightforward. 

1. Introduction 

1.1. The present paper is concerned with the study of the vibration of a thin elastic 
plate with free edges. Let us consider a plate that occupies in its natural condition a plane 
domain D in the xy-vl&ne bounded by a sufficiently regular contour C. It is known [1,2] 3 
that the problem of determining the frequencies of this plate may be reduced to the following 
eigenvalue problem, denoted by (Pr.), for the differential system consisting of the differential 
equation 



A 2 w 



/d 2 . d 2 \ 2 



-\w 



hxD, 



(1) 



with the boundary conditions 



dAw . d d 2 w n 

-^— + (1— \x)-r v^7 =0 on C> 
dn as dndt 

nAw+ (1 -M) dn2 =0 on C, 



(2) 



where /i is an elastic constant, called Poisson's ratio, such that 0^ju<3 (the mathematical 
theory applies for ^m<C1)> & n d where s is the arc length parameter of C and b/dn, b/dt repre- 
sent respectively the derivatives in the directions of the outer normal vector n and the tangential 
vector t to C at the point under consideration. In the case of angular points on C, where the 
direction of n varies discontinuously with a jump, w is subjected to an additional boundary 
condition [1]: b 2 w/dndt at a point P on C tends to the same value as P approaches any one of 
the angular points from either side. 

The main object is to propose a variational method that enables us to calculate approxi- 
mate values of the eigenvalues with a rigorous estimate of the error; in other words, to calculate 
upper and lower bounds for the eigenvalues. Such a method seems interesting and important 
from the viewpoint of application, inasmuch as it appears, even for most elementary shapes of 
C, almost impossible to obtain an explicit solution of (Pr.). On the other hand, usual approxi- 
mating methods applied to (Pr.) seem to suffer from the lack of error estimation. For example, 
the well-known method of Ritz is applicable and was actually applied to the case of a square 
plate [2]. Although we know empirically that his method provides us with approximate 
values whose accuracy is satisfactory in most cases dealt with in practice, it should be noted 
that they are only known to be upper bounds for the quantities in question and nothing can 
be said about the error bounds. See [10] for a comparative discussion of existing methods, 
some of which furnish lower bounds. 



i University of Tokyo, American University, and National Bureau of Standards. 

2 Department of Physics, University of Tokyo. 

3 Figures in brackets indicate the literature references at the end of this paper. 

169 



1.2. Our method is based on the following theorem, proved in [3], concerning a self- 
adjoint operator in a Hilbert space 4 of the type T* T, which reads, in a form slightly modified 
for our later convenience, 

Theorem l. 5 Let T be a closed linear operator with its domain dense in a Hilbert space 
<Q and its range in a second Hilbert space Q' and let T* be the adjoint of T. Let a 2 <C\<ifi 2 
(0^a</3) be an open interval containing at most one nondegenerate eigenvalue oj the operator 
T*T but no other points oj its spectrum. If u and v are elements belonging to the domains of T 
and T*, respectively, such that 

\\u\\t*0, !M|t*0, (Tu,v) = (u,T*v)^0 (3) 

and 

(v-*)(P-v)>e 2 (4) 



hold, where 

_ (Tu,v) 



\u> • \\v\\ 



CWTuW 2 II T*v\\ 2 1 
\\iu\\ \\i v\\ 2 

\\u\\ 2 + HI 2 J * ' (5) 



then there exists certainly an eigenvalue X' of the operator T*Tin the interval (a 2 ,P 2 ) and the following 
inequality is valid: 6J 

L^X'^U, (6) 

where 

£=^)=(*-~) 2 , U=U(u,v)=(r,+^y. (7) 

In order to make this theorem applicable to (Pr.), it is required to construct the spaces 
£>, §' and define the operators T,T* in such a manner that the eigenvalue problem T*Tw=\w 
is equivalent to (Pr.) and the quantities in (5) are practically calculable, at least for the u and 
v chosen. (Pr.) is nothing but the eigenvalue problem for a self-adjoint operator H (in a real 
Hilbert space § of square integrable functions on D), defined as 

Hw=A 2 w (8) 

for w subjected to the same boundaiy conditions as those of (Pr.), in addition to certain regu- 
larity conditions making A 2 w well-defined. Therefore, our problem is to reduce this pre- 
assigned operator H to the form 

H=T*T (9) 

with a suitable choice of § r , T, and T*. Here, u and v may be called "trial functions," or 
possibly " trial vector functions," in estimating X', a better choice of which will yield a better 
estimate. Finally, a 2 and fi 2 are, loosely speaking, a rough upper bound of the next lower 
eigenvalue to X' and a rough lower bound of the next upper eigenvalue, respectively. Such 
rough estimates will be obtained by other methods (for instance, the Ritz method for a 2 ), or 
occasionally by solving an appropriate auxiliary problem, as in the case for /3 2 concerning a 
square plate, treated later. The decomposition of H into the form T*T is carried out in 
section 2 for a general shape of C. This is followed by some supplementary remarks on theorem 



4 Concerning the theory of Hilbert space, we mostly follow the notations and terminology of Stone [4]. For operators between two Hilbert 
spaces, see Murray [5]. 

8 An elementary proof of this theorem will be given elsewhere. 

6 Any eigenvalue of T*T is nonnegative, because (T*Tw,w) = \\ Tw|| 2 ^0. 

7 v — — >0 and 7/H >0 follow from (4). 
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1 for improvement of the approximation. In section 3, confining ourselves to the case of the 
square plate bounded by x= ± 1, y= ± 1, we illustrate in detail the method of actual calculation 
along the lines mentioned above and give a numerical estimate of the smallest positive eigen- 
value of (Pr.), taking ju=0.225, with a relative error bound less than 1/1000. Section 4 is 
devoted to concluding remarks, with a brief reference to other methods with the same object 
that are applicable to (Pr.). 

2. Decomposition of the Operator and Remarks on Theorem 1 

^ The first part of this section is devoted to decomposing H, the differential operator of (Pr.) 
including the boundary conditions, into the form T* T in the sense stated in the introduction, 
namely, to constructing T and T* so that H=T*T, i. e., T*Tw=\w, implies (Pr.). Strictly 
speaking, the differential operators appearing here as well as in the sequel must be replaced by 
their closed extensions with certain extended domains, or in other words, differentiations should 
be interpreted in a certain generalized sense. Nevertheless, detailed discussions of this situa- 
tion will not be given, for they seem to be unnecessary for practical application, which is the 
ultimate object. 

2.1. We begin with the following heuristic considerations. 8 Let T be a closed operator 
from § into €>' and T* its adjoint. Then the eigenvalue problem T*Tw=\w is equivalent . al 
least under a suitable assumption on the spectrum of T*T, to the variational problem 

8\\Tu\\ 2 =o, ||^|| 2 =i ? 

where 6 is the first variation. Thus, T*Tis the " gradient operator" of the quadratic functional 
||7V|! 2 in the sense that 8\\ Tu\\ 2 =2(T*Tu,bu). On the other hand, as explained in textbooks 
[1, 8], (Pr.) is derived from the variational problem 



5J[u] = 0, I uHxdy = 1 , 

-J n 



where J[u], which represents (apart from a constant factor) the elastic energy of the plate, is 
expressed by 



^#«--) {&$ -(&)'}]**. 



(10) 



that is, II is the gradient operator of J[u]. Accordingly, § being the same as that introduced 
in defining 77, it is reasonable to adopt $' and Tsuch that 

\\Tu\\ 2 = J[u] (11) 

holds. In fact, we proceed as follows. Let L 2 (D;p) denote the real Hilbert space consisting 
of functions integrable (in square) on D with norm defined by 

|k|| 2 = pj u 2 dxdy, 

where p is a positive constant. 

In view of (11) and the fact that J[u] can be written as 

J M ^£[MA^+(l-,){(g)V(g) 2 +2( d ^ / ) 2 }]^ (12) 

we define §' by 

i3' = J L-(^;M)Xi 2 (/>;l-M)Xi> 2 (L';l-M)Xi 2 (L»;2-2 M ), (13) 



8 Some general devices useful in such decomposition are proposed in [6]. 
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where X represents the Cartesian product of Hilbert spaces. It should be noticed that an 
element v of §' is a four-component vector function v={vi,v 2 ,VzjV4\ defined on D. Now the 
requirement J[u]= \\Tu\\ 2 is evidently fulfilled by the operator T, defined by 

Tu=Tu(x,y)={Au, £* 0, *L} &'. (14) 

Elements u belonging to the domain 3) of T are subjected to certain regularity conditions, but 
are otherwise arbitrary, i. e., free from any boundary conditions. Though these regularity 
conditions are too complicated to be specified explicitly in elementary terminology, it will be 

A 

enough for our purposes to know that the class C of functions u, such that u, bu/bx, bu/by are 
continuous and b 2 u/bx 2 f b 2 u/bxby } b 2 ujby 2 are piecewise continuous on D+C, is contained in 

S); namely, &($). 

T* is defined (together with its domain 3)*Gp') by the condition that 

(Tu,v) = (u,T*v) (15) 

be valid for any w££). To determine the form of T 7 *, we note that 

(Tu,v)=jJ Lli (Au)v 1 +(l- f i)^v 2 +(l-„) |^3+2(1 - M ) ^ njixdy. (16) 

In reducing the right side of (16) to the form u-T*v dxdy, we can resort to devices similar to 

those usually employed in deriving the differential equation and the boundary conditions of 
(Pr.) from the variational problem concerning J[u], Before doing this, we introduce some 
notations, some of which appeared in section 1. 
Let 

s= arc length parameter, 

n=n{s) = (n z ,n y ) =unit outer normal vector to C, 
t=t(s) = (t x ,t y )=umt tangential vector to G in the positive direction, 

and note that 

d 5 , b b h , b 

K-=n x -zr — \-n v zr-y "=zi=tx si — Kv s^" , (1«) 

bn bx v by bt bx y by 

and also 

bx= n *Yn +t *bt' Vy^Yn^bt (18) 

Hereafter it is assumed that the components of n and t are differentiable except at possible 
angular points Pi,Pz, . . -,P n with s=s 1 ,s 2 , . . .,s n . Also it is known that, for any suffi- 
ciently regular functions F(x,y), G(x,y), the identities 



f |j^ 6dxdy=- f F ^ dxdy+ f F-Gn x ds, 
f ^ Gdxdy=- \ F ^ dxdy+ f F-Qn y ds, 



(19) 



are true. 

Provided the components of v are sufficiently regular that all procedures used below may be 
justified, the expression for (Tu,v) in (16) becomes, by some small calculations that make use of 
(18) and (19), 
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where 



(^)=X-[mA, 1+ (1-m){0+^+2 ^yy+-Ju**+f]» Qds+J^Bds, 

}1 



P=P(V,8)=^+ 



%«*%»*%■■ 



Q=Q (p,8) =iaVi + (1 — n) { v 2 n 2 x +v 3 nl + 2vji x n y } , 

R = R(V,8) = (1— M) {»2Wx<z + tl3WA + »4(^+V*) 1 

Furthermore, because at a point P with s=<t we have 



/cto\ fdu(8)\ 
\i>tJ P \ ds ) 9m9 



and consequently 



^rr Rds= -j- jBrfs, 



we can rewrite the last term of (20) by partial integration, taking account of the fact that C is 
closed and that R, as well as the components of n and t, is differentiable with respect to s except 
at s=Si, 



., s n , so that 



(22) 



where D*=lim R—lim R, 



T-T Rds=— u— ds—^u(si)D i} 

JcOl Jc (IS i = i 

(1=1,2, . . .,n). 
The substitution of (22) into (20) results in 

(23) 

This expression should be equal to (u,T*v). Because u and bu/dn on G are arbitrary, we 
immediately conclude that this is true if and only if v satisfies the boundary conditions 

d 



P +ds R =°> 
Q=0, 

D t =0, (i=l,2, . . . , n). 



(24) 



Then T 7 * is given by 



^r* {W8) ,}=M,+(i-,){^+|+2^}. 



(25) 



Though an element v belonging to 55* must obey certain regularity conditions besides the 

boundary conditions (24), we will not discuss these in detail. Again it is sufficient to know 

A 
that any v with components belonging to C and satisfying the boundary conditions (24) is 

contained in S*. 

The construction given above shows that the operators T and T 7 * are certainly the desired 

ones, but it may be worthwhile to verify directly the equality H=T*T, i. e., the equivalence 

of the eigenvalue problem T*Tw=\w to (Pr.). 
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Let w belong to the domain of T*T. Then Tw belongs to J)* and consequently is subjected 
to (24). But 



dAw fdh 

" (1 ~ M) \dr 



w 



dn 

dAw 
dn 



. d z w 
n x +^-^n y + 



d d w 



dx s ' l ' dy s ' vy ' dx 2 dy ' 



d 3 w 
dxdy 2 



■} 






?)"•}] 



dAw , /i .5 /d 2 w . d 2 w\ 



= ^ Aw 

on 



by virtue of (17). Similarly, we obtain 



<2(Tu>)=mAw+(1-/«H ^r~ 



'{ 



d 2 w 



n x + ^TJ n l + 2 



d 2 w 



dy 2 v dxdy 



n x n 



=fiAw+(l — ij) 



d 2 w 



d 2 w 



<{ 



dn 2 

d 2 w , , d 2 w 



fluty 



dx 2 xx ' dy 2 vv 



=(1-M) 



d<dn 



} 



d 2 w 
dxdy 



(n x t y +n y t x ) 



} 



Thus, w satisfies all boundary conditions of H and therefore belongs to the domain of H. Con- 
versely, by a similar argument it is seen that any element in the domain of H belongs to the 
domain of T*T. Finally, we observe that 



dxdy \dxdyj J 



= A 2 w, 



as it should be. 

2.2. The following remarks on theorem 1 include some devices that would be of use in 
obtaining a good estimate, particularly in constructing a desirable pair u and v. 

Remark 2.1. The estimates given by theorem 1 are sharp when e 2 is small. But it is 
known [3] that e 2 can be written in an alternative form, 



\\u\\ 



\W 



+ 



and e 2 = occurs if, and only if, 



1U 11*11*' 



and 



\u\\ 



1 v M\* 



This implies that T*Tu=rfu, so that rj 2 is equal to an exact eigenvalue X' of T*T and u is the 
corresponding exact eigenfunction w. Thus we may expect a good result when u and v approxi- 
mate, apart from constant factors, w and Tw closely. This interpretation of u and v should be 
taken into account, in particular when some knowledge concerning the behavior of w is given 
in some way or other, for instance from physical considerations or experimental results. As 
for precise discussions of the relations between u, v, and w, reference is made to the paper by 
Kato [3]. 
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Remark 2.2. If u and v are not completely preassigned but admit of some free parameters, 
it would be intrinsically the best way to choose those parameters first to minimize f/and then 
maximize L. In practice, however, it appears sufficient as well as convenient for calculation to 
determine them by minimizing e 2 . A standard form for such adjustable trial functions is 



U= aiU a) + a 2 U (2) + . . . +a m u (m) , 
v=/3 1 v^+(3 2 v™ + • • • + ft^ (w) , 



(26) 



where u {1) , u (2) , . . .,u {m) are fixed functions belonging to J) and v {1 \v i2 \ . . .,v (n) are fixed vector 
functions belonging to 55 *, and where a h a 2) . . .,a m , ft, ft, . . -,ft are Iree parameters. If we 
calculate the quantities in (5) with these u and v, then \\u\\ 2 , ||7^|| 2 become quadratic forms in 
{a z ), \\v\\ 2 , \\T*v\\ 2 quadratic forms in {ft}, and (Tu,v) a bilinear form in {a { } and {ft}. Thus e 2 
becomes an algebraic expression in those parameters whose minimum is desired. There is 
little difficulty in the necessary computations, especially when we can resort to an automatic 
computer, even if m and n in (26) are not very small. 

Remark 2.3. Suppose that a closed subspace 2ft of £> reduces the operator H=T*T, 
that is, let HweWl for any w belonging to the intersection of Tl and the domain of //. Then we 
can regard H as an operator in the Hilbert space $Jl denoted by H m . It is known that eigen- 
values and eigenvectors of H^ t are also eigenvalues and eigenvectors of H, so that we may speak 
of an eigenvalue of //in 9JJ. If the eigenvalue X' of //in question turns out to be an eigenvalue 
of H% h we can derive an estimate for X' from theorem 1 more conveniently by restricting T to 
9K, and T 7 * to a suitable closed subspace Wl* of ,£/, with the properties 

218Hn£)c2R*, r*(SW*n©*)c3W. (27) 

This procedure is particularly advantageous when a degenerate eigenvalue of the H in question 
is converted into a nondegenerate one of Z/g# by a suitable choice of 3JJ so that applica- 
tion of theorem 1 becomes possible, or when the neighboring eigenvalues of H m pre more widely 
separated than those of //, thus permitting us to take smaller a and larger f$ so as to make the 
estimate better. Moreover, though the smallest eigenvalue of H in our problem is equal to 
zero, the associated eigenf unctions being of the type ax-{-by-\-c, we can sometimes find 9ft such 
that the eigenvalue X r , originally intermediate, becomes the smallest eigenvalue of II^\ and 
hence its upper bounds can be obtained by 



tf,«HffcV f fceaKn©), 



(28) 



the so-called Rayleigh principle. This seems more convenient than the estimation by the right 
inequality of (6), for it does not involve v } an element of J)*, for which a good choice is not easy 
because of the complication introduced by the boundary conditions. Therefore, we may 
prefer U R to U as an upper bound for such an eigenvalue. 

Nontrivial 9K, 3JJ* as mentioned above are often obtained by considering the symmetry 
properties of the operators H, T, and J 7 * reflecting those of C, as we shall illustrate in the 
succeeding section in the case of a square plate. 

3. Vibration of a Square Plate 

In this section the general results obtained so far are applied to the case of a square plate. 
Let C be the square with vertices (1,1), (— 1, l), (—1,-1), and (1,-1), which is the same as 
that treated by Ritz [2]. Bounds will be calculated for the smallest positive eigenvalue of H= T*T 
according to theorem 1 and the Remarks in section 2, i. e., lower bounds by means of (6) and 
upper bounds by means of (6) or (28). 
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3.1. We begin with a preliminary consideration. Let 

J[u]=(l-ii)\Jul x +u$ y )(fady. 



(29) 



Let the self-adjoint operator H in § be the gradient operator of J[u]. Explicitly, anj r element 
w in the domain of H is subjected to the boundary conditions 



d 2 w_d 3 w_Q 
da; 2 da; 3 



(x=±l,-l<y<l) 



^=^=0, (-l<*<l,y=±l), 
dy 2 dy 3 " 



(30) 



and we observe that 

Hw = (1 ~4^ + W (31) 

The eigenvalue problem H w=\w is reduced, by separation of variables to the eigenvalue prob- 
lem 



d A u 



(32) 



u"(±l)=u"'(±l)=0, 

and a similar one in the variable y. Equation (32) is nothing but the eigenvalue problem treated 
in the analysis of vibrations of a, free bar of length 2, and is completely solved in textbooks 
[7,8]. If we denote the eigenvalues by k* (n=0,l,2, . . .) and note that is the unique degen- 

nun t n y-vr-» i-\ ttti + In r\i m»/~vY-» tit v-» /-»+■ i /-»v-» c-t n r\np tita Kotta 



erate one with eigenf unctions a + bx, we have 

h=h=o, o<k 2 <k 3 <:k 4 < . . 

where 9 

& 2 =2.3650204, fc 3 =3. 9266023, 
& 4 =5.4978039, k 5 = 7.0685828, 

The normalized eigenf unction u n =u n (x) associated with ki is given by 
1 



(33) 



Ua=-=> 



V2 



U ^h 



u n = 



cosh k n 'Gos Jc n x-\- cos Zvcosh k n x 
Vcosh 2 k n +cos 2 k n 

sinh k n • sin k n x+$m k n • sinh £ n £ 



Vsinh 2 k n — sin 2 & ra 



for even w>2, r 



for odd n > 3. 



(34) 



It should be noticed that u w is an even or odd function, according as n is even or odd. 
The eigenvalues and eigenfunctions of H are given by 

Xm»=(l — n)(ki+ki), w mn =u m {x)u n (y), (m,n=0,l,2, . . .). (35) 

We note that w mn (m, 71=0,1^2- . . .) constitutes a complete orthonormal system of eigenfunctions 
of the self -ad joint operator H. 

3.2. Comparing two quadratic functionals J[u] and J[u], or two operators H and H, 
we can not only determine the constant p in theorem 1 in a rigorous fashion, but also the sym- 

9 fc n is obtained as a root of equation tan fc=tanh k or tan k= — tanh k, according as n is odd or even. 
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metry of the eigenfunction associated with the smallest positive eigenvalue. We note that 

J[u]^J[u]. (36) 

From (36) it follows, by virtue of the properties of H stated above, that the sped nun of H 
consists only of discrete eigenvalues, and also that, if we number all eigenvalues of H and II in 
ascending order, counting any degenerate one repeatedly by its multiplicity, we have 



A N ^A Nf 



(37) 



where A^ and A^ are the iVthjeigen value of H and the iVth eigenvalue of H, respectively [7]. 
Thus the eigenvalues (35) of H are the lower bounds of the corresponding eigenvalues of H, 
and hence we can take (3 as 



:(1 -„)(% + «), 



(38) 



if the eigenvalue X' of H in question is the rth eigenvalue from below and \ mn is the (r+l)th 
eigenvalue of H. 

We turn to the study of the symmetry properties of H and H, i. e., of closed subspaces 
specified by the symmetry properties of their elements, and to the reduction of each of these 
operators. For the square plate under consideration, the boundary conditions of // follow 
from (2) and the remark there and take the form 



W+'^wh' 






d 2 w . d 2 w 



d 2 w . d 2 w 



d 2 w 



dxdy 



=0, 



(x=±l,-l< y <l), 

(y=±l,-l<x<\), 

(x=±h-\<y<i), y 

(y=±l,-l<x<l), 
(s=±l,y=±l). 



(39) 



In the deduction of the last condition of (39), use is made of the fact that, for instance, 

d 2 w d 2 w ,, . , H d 2 w d 2 w 



dndy dxdy 
and hence that the condition 



,. d 2 cn ,. d 2 w 
lim ^ ^, =Jim: 



on the side x=l, 



' dndt 



dndt 



means 



dndt 



d 2 w 



dxdy 



on the side y=l, 



z = l 



x-*l 
y=\ 



dxdy 



=0 at the vertex (1, 1). 



Owing to the symmetry properties of (39) and those of the differential operator A 2 , H is 
reduced by the following subspaces: 

Wl(o,o) = set of functions odd both in x and in y, 

Wl(o,e) = set of functions odd in x and even in y, 

<jffl(e,o) = set of functions even in x and odd in y, 

dlt{e,e) = set of functions even both in x and in y, 
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(40) 



which are mutually orthogonal and whose direct sum is the whole space fQ. We see immediately 
that iJis also reduced by each of these subspaces, so that the inequality (37) still remains valid, 
if we regard H and H as operators in any of them. Let the eigenvalues of H in these subspaces 
be denoted as 10 

\><Xi<X,<. . . , in2W(o,o), 

Mo<Mi<M2<- • • , either in $l(o,e) or in $Jl(e,o), 

*>oOi02<. . • , in Wl{e,e) 

the union of which is the total spectrum of H. On the other hand, the eigenfunction w mn of H 
with the eigenvalue \ mn belongs to 9J?(o,o), $Jl(o,e), 2JJ(e,o), or SBl(e,e), according as (m = odd, 
w=odd), (ra=odd,7i=even), (m=even, n=odd), or (m = even, 7i=even). 

Because the zero point set 9t of H (the eigenspace with the eigenvalue 0) consists of func- 
tions of the type ax + by+c, it follows that 

m(o,o) a.% 

9W(o,e)^9?={ u; u=ax} 
yjl(e,o)^yi={u;u=by} 

and hence that 

In consequence, the smallest positive eigenvalue of H is equal to min(X ,Ad^i). We can 
determine the order relations between X , mi, an d ^i as follows. Because (40) means that /x 
is degenerate neither in 3R(o,e) nor in $Jl(e,o), and v is not degenerate in d)l(e,e), we note, 
according to (37), that \i2( = \2i) and X22 are lower bounds of ^i and v Y , respectively, and hence 
we have, by means of (33) and (35), 

Mi^Xi2 = X 2 i=(1-m)£ 4 2=(1-m)X31.28 . . . ,) 

\ (41) 

^\22=2(1-m)£!=(1-m)X62.59. ... ) 

On the other hand, the substitution of the trial function u = xy^Wfl(o,o) into (28) yields 

X ^(1-m)X18. (42) 

Combining (41) and (42), we conclude that the smallest positive eigenvalue of H is the smallest 
eigenvalue X in 9K(o,o). Moreover, owing to (42) and the relations X n =0, X 13 =X 31 = (1 — /*)X 
237.72 . . ., we obtain 

Xii<Xo<X 13 =X 31 , (43) 

which shows that X n is the unique eigenvalue of H less than X , and hence X is not degenerate- 
Also we have, by means of (37), X 13 ^ \ h which enables us to take as $ in theorem 1, with respect 
to X' = X , any number such that 

2 ^(l-iu)X237.72 .... (44) 

Furthermore, we divide Wfl(o,o) into the following two mutually orthogonal subspaces 

m and ar, 

W={u; u £ 2W(o,o), u(x,y) =u(y t x) } , 
W = {u; u £ 2W(o,o), u(x,y) s —u(y,x) } , 



10 Evidently, the spectra of Htyfl ( ,e) and those of Hffll («,„) are identical, and any eigenfunction of one of them is converted into an eigenfunction 
of the other by permutation of i and y. 
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each of which is seen to reduce both of H and H. 

From (43) it follows that X is an eigenvalue of II not in 2ft' but in 3ft, because the smallest 
eigenvalue of 27 in 9JT is X 13 with eigenfunctions Const- (w i3 — W31), though (his is at the same time 
an eigenvalue of // in 5Df with eigenfunctions Const -(wiz+ w 31 ). In calculating hounds of X , 
to which we confine ourselves hereafter, we may therefore restrict H and T to 2ft. 

3.3 We go on to the construction of the trial elements n u and v in theorem 1. Introducing 
six elements u a) , . . . ,u {6) of 2ft, given by 

tt a, =Pi(*)Pi(»), u*>=P l (z)P t (3/)+P t (x)P l (y), 

u™=P z (x)P z (y), ^^PiW^W+AWAW, 

u™=P z (x)P 5 (y) +P 5 (x)P z (y), u™=P 5 (x)P 5 (y), 

where P t represents Legendre's polynomial of degree i (i= 1,3,5), we put, in accordance with 

(26), 

u=a l u {1) + a 2 u i2) + . . . +cx 6 ?/ (6) (45) 

with free parameters «i,a 2 , . . . ,« 6 to be determined later. This u evidently belongs to 
3B n J), because it is a polynomial and surely satisfies the regularity conditions of 2). 

The other element v= {v u v 2 ,v z ,v 4 } is constructed as follows. We can easily verify that the 
boundary conditions (24) of £)* now reduce to the following forms: 



^ +(1 " m) ^ +2(1 - m W = ' 

/Xt'i+(1— /X)« ? 2 = 0, 
^i+(1— M>'3 = 0, 



(x=±l, -l<y<l), 

(-l«lj=±l), 

(z=±l, — l<y<l), 
(-1<«<1, 2/=±l), 
(*=±1,»=±1). 



(46) 



In view of the symmetry of 2ft, on the other hand, we can take as 2ft* the subspace com- 
posed of elements r with the following properties of symmetry: 



Vi, v 2 , v z are odd both in x and in //, 

v 4 is even both in x and in y, 

v x (x,y) =Vi(y % x), v 2 (x,y) =v z (y,x), v A {x,y) = v 4 (y } x). t 



(47) 



Then the relations (27) obviously hold. 

A general element v satisfying the conditions (47), and whose components are polynomials 
in x,y with the degrees shown below, is expressible as 

0={t>l,«W4}, 



t>i = a x xy + a 2 (xy z + x*y) + a 3 ;cV + a 4 (x?/ 5 + x 5 y) , 
2? 2 = b x xy J r boxy* + ft^fy + &3#V + b *xy 5 + 6 4 ^ 5 2/, 
v z =b x xy + 6^2/ 3 + &2« 3 2/ + b z x*y z + 6 4 ^2/ 5 + b 4 x 5 y, 
v 4 =Ci + c 2 (x 2J r y 2 )+c< i x 2 y 2 +c i (x*+y*). 



(48) 



11 If we are contented with a rough estimate of X°, we can construct the trial elements more easily (see [9]). 
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Moreover, for such a trial function v the boundary conditions (46) are equivalent to 

bv 1 bv 2 dv A 



a ^ — hs 
dr 5/ 



where <r = /z/(l— /*). 



dy 

av ] +v 2 = 0, 
v 4 =0, 



(x=l, -l<y<l), 

(3=1, -1<»<1), 



(49) 



^(«l+«2 + « 4 ) 


+(61+62+^) 




=0, 


a(a 2 + 3a 3 ) 


+ (6 2 +36 3 ) 


+ 8c 4 


=0, 


^(«2 + « 3 ) 


+ (62+63) 




=0, 


ca 4 


+64 




=0, 



Substituting (48) into (49), putting x=l and equating the coefficients of y, y z , y 5 to zero 
or putting x=y=\, we obtain the conditions for v given by (48) to satisfy (49) and hence to 
belong to £)*. The results are the following six equalities: 12 

<T(a 1 + 3a 2 + 5a 4 )+(6 1 + 362 + 56;) + (4c 2 +4c 3 ) = (0 



(50) 



c 1 + 2c 2 +c 3 +2c 4 =0. 

By inspection, we find the following 8 independent families, (1), (2), . . . , (8), of 
parameters satisfying the system (50) of 6 equations in 14 unknowns: 



d=l, C2=-l/2, 

£l~ ^—l? ^2 = 1? 

c, = l, c 3 = — 1, 

d = l, c 4 =-l/2. 

The absent members of each family are equal to 0. 

In this way, we have the following 8 elements v (Jc) ={v[ k) , v { 2 \ v^, v {k) \ (Jc=l, 2, . . . , 8) 
contained in 9JJ*^S*. 



(1) 


«i=l, 


6,= — (7, 


(2) 


02=1, 


62=62= — (T, 


(3) 




6,= -1, 6 2 =1, 


(4) 


a 3 =l, 


6 3 = — <r 


(5) 






(6) 


a 4 = l, 


64=64= — cr, 


(7) 




6 1= -1, 64=1, 


(8) 




6 2 =-2, 63=2, 



: {sy, 



-crary, 



- <rai/, } , 



v (2) ={xy 3J rx*y, — a(xy ZJ rx 3 y), —<r(xy 6 -\-z*y), 0} , 
z> (3) = {0, —xy+x z y, —xy+xy 3 , l-(x 2 +y 2 )l2}, 
v ( * } ={ x*y* , — ax*y 3 , — ax 3 y* , } , 
^> = {0, 0, 0, (l-x 2 )(l-y 2 )}, 
0< 6 >={ay 5 +arty, — <Ka# 5 +s 5 2/), —<r(xy 5 +x 5 y), 0} , 

^ (7) =-{0, —xy-^x^y^ —xy + Xy 5 , 1 — £ V } , 



,(8). 



;0, 2(-a7y 3 +aty 3 ), 2(-rV+rV 3 ), l-(^+?/ 4 )/2}. 



12 The result of equating to zero the coefficients of y 5 in the first and second conditions of (49) coincides and gives the fifth condition of (50). 
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With these 8 elements as basis, we construct the trial function r in accordance with (26) : 

»=ft» (1) +A»®+ • • • +&*' (8) (51) 

with free parameters &i,p 2 , . . . ,/3 8 . 

Making use of the trial functions u and v given by (45) and (51) respectively, we get, after 
the necessary integrations, 

6 

II^H 2 =Z) Ai^a^il-^^A^a^ (A ik =A ki ), 

i,k=l 

6 

IMF=S BuPLipi k , (B ik =B ki ), 



t*v\\*=j: ftAft=d-M)2:a'A&, (c«=cw, 



».*-i 



IHI*=S #«ft&=(l-/i)X^;*P*&, (Da=D kt ), 



t = l fc = l 



with i4«, 5 iA; , CI*, D' ik , Ea given by the tables. 





1 


2 


3 


4 


5 


6 


1 


8 


16 


8 


16 


16 


8 


2 




400(7 + 312 


96 


L120<7+672 


192 


96 


3 






600 , 2616 

7 7 


360a + 96 


240(7+816 


288 


4 








4648o- + 3336 


1008(7 + 336 


240 


5 










22248 . 37296 
11 a+ 11 


1680 , 17520 
11 a+ 11 


6 












840(7+2640 



#;/ 



l 


1 


2 


3 


4 


5 


6 


B<*=0 


(*VJfe) 


B« 


4 


8 


4 


8 


8 


4 


9 


21 


49 


33 


77 


121 
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CL 



i 


3 


3 


3 


3 


5 


5 


5 


7 


7 


8 


For other values 
of i and k, 

c;*=o. 


k 


3 


5 


7 


8 


5 


7 


8 


7 


8 


8 


W/fc 


64 


178 
3 


256 
3 


348 
5 


256 
9 


512 
9 


256 
5 


8704 
63 


4096 
35 


17664 
175 



D[ k 



i \ 

\ 


1 


2 


3 


4 


5 


6 


7 


8 


1 


4 
9 r 


8 
15 T 


16 
45* 


4 
25 r 





8 
21 r 


32 
63* 


32 

75* 


2 




368 
525 r 


64 
175" 


8 
35 T 





496 
945 T 


512 
945* 


256 
525* 


3 






144 
35 


16 

175 


128 
45 


32 
135* 


992 
189 


2528 
"525 


4 








4 
49 r 





8 
45 r 


32 

225* 


32 

245* 


5 










512 
225 





256 

75 


1024 
"315 


6 












656 
1617 7 


192 
53^* 


64 
189* 


7 














119744 
17325 


1088 
175 


8 
















63424 
11025 



E ik 



i 


1 


1 


1 


1 


2 


2 


For other values of i and k, 
E' ik = 0. 


k 


3 


5 


7 


8 


7 


8 


E ik 


16 
3 


32 

9 


64 
9 


32 

5 


64 
21 


64 
35 



Here a and r are constants given by <r= 



l-V 



1 + /X 

"(I-/*) 



In the numerical computation we take ^ = 0.225, the value used by Ritz. Then the quan- 
tity 2e 2 becomes an algebraic expression of the form 



_(x^Ey) 2 _ 
x'Bx^y'Dy (x'Bx) ( KV f Bf! 



2 _x'Ax , y'Cy 
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(52) 



where x, y are the column vectors with components [on], {&}, respectively, x', y' are their 
transposed vectors, and A, B, C, D, E are matrices with numerical elements. 

In order to minimize e 2 , we choose the following procedure. First we minimize the first 
term on the right-hand side of (52). This gives rise to an eigenvalue problem of order 6, for 
which we determine the lowest eigenvalue m (which is the mini mum value in question) and 
the associated eigenvector x=x u Then we put x=Xi into (52); this converts (52) into a frac- 
tional form of which the numerator and denominator are known quadratic forms in y. We 
minimize this expression; this is equivalent to solving another eigenvalue problem of order S 
and we determine the lowest eigenvalue m' and the associated eigenvector y = y\. In the next 
step we should put y = yi into (52), converting (52) into a ratio of two quadratic forms in x, 
and determine its minimum m". This procedure could be repeated indefinitely, and the 
sequence m' ', m" , . . . is certainly nonincreasing. It is not clear whether these values con- 
verge to the true minimum of 2e 2 . But this is not important, for actually we need not obtain 
this true minimum; it is sufficient for our purpose that a very small value of 2c is given by 
some set of x and ?/. 

The actual computation was done by using SEAC, and it turned out that the value of 
m r was sufficiently small, so that further computation was not necessary. To avoid the uncon- 
trollable error arising from the use of the machine, the components of ./',,//, thus obtained have 
been rounded off to five significant figures, and the value of 2e 2 for these arguments has been 
calculated anew on a desk calculator. The final values arc 4 

«!= 3.6100 000, ft = 2.0362 000, 

a,= -0.2071 000, 2 =- 1.1853 000, 
0:3=— 0.0357 000, 08= — 1.6269 000, 
a 4 = 0.0335 000, 4 = 0.3121 000, 
a 5 = -0.0012 000, 5 = 0.5907 000, 

a 6 =— 0.0017 000, 6 = 0.0112 000, 

7 = 0.0176 000, 

8 = 0.1233 000, 

e 2 = 0.009172, rj= 3.6266247. 

We can take a=0 and 0=13.5 in (7) for U and L in conformity with (44). Thus we 
obtain 

£7= 12.455433, L= 12.430613, 

where the rounding-off errors are within 5X10 -5 by the most conservative estimates. 'Finis 

12.4305<X <12.4555. 

For the upper bound we can also use U R , given by (28). This gives 

£7^=12.454620, 

which is somewhat sharper than U. In this way we arrive at the final result : 

12.4306<X <12.4547. (53) 

If we adopt X = 12.4427, the average of upper and lower bounds, as the smallest positive 
eigenvalue of H=T*T, the relative error is less than 1/1000. 

The corresponding value of the fundamental frequency of a thin elastic square plate 
(jj. = 0.225) with free edges is 



^=J-\=0.3420xJ^\, 
l-fx 2 V P a \ P a 2 



2V3ttV 

with relative error less than 1/2000. Here 2a = edge length, 2h= thickness, p = density, and E= 
Young's modulus. 
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A remark is necessary regarding the relation between the result of Ritz [2] and that of this 
study. Ritz gives the estimate X ^ 12.43, and this appeared exceedingly sharp in view of the 
results shown here (53). Ritz used products of eigenf unctions of free bars as trial functions 
u {i) . This suggested that these trial functions might lead to a sharper upper bound than the 
one presented in this paper. Because the last figure in Ritz' result is unreliable, the same 
calculation was carried out for purposes of comparison by retaining more significant figures. 
The result obtained gives only A <12.488, a less sharp estimate than the one herein. This 
seems to show that, as far as the fundamental frequency of a square plate is concerned, use of 
polynomial trial functions is preferable to products of eigenf unctions of free bars. 

4. Concluding Remarks 

4.1. This method is also applicable to a plate whose edge is not wholly free but subjected 
to conditions of a more general type. For instance, let us assume that C> the boundary of the 
plate under consideration, consists of three arcs, C , C h and C 2 , and that the edge is free, sup- 
ported, and clamped along C , C h and C 2 , respective]y. Among all angular points on C, we 
denote those on C by Pi, P 2 , . . ., P n with s—Si,s 2j . . .,s n . Then as in the case of the free plate 

discussed above, the vibration of this plate is characterized by the eigenvalue problem of a 

^ . a A 

self-adjoint operator H in §. Again, H is the gradient operator of a quadratic functional J 

[u]j which is the restriction of J[u] to the class of functions u satisfying the boundary conditions 



^=0 on Ci + (V 






_ =o on C 2 . 



(54) 



A A A 

Thus, following a course parallel to that of section 2, we can reduce H to T*rand make theorem 
1 applicable. The results of this decomposition are as follows. § and §' being the same as 

before, Tu and T*v are formally identical with Tu and T*v, respectively, and now any u be- 

A 
longing to the domain of T is subjected to the boundary conditions (54) and any v belonging to 

the domain of T* is subjected to the boundary conditions 

Q=0 on Co, +C l9 

P+jR=0 on Co, 

D t =0, (1=1,2, . . .,n), 

where P, Q } R, D t are those defined in section 2. 

Applications of theorem 1 to the cases of such mixed boundary conditions will be treated 
elsewhere. 

In addition, if both C and Ci are empty, i. e., the whole edge is clamped, another de- 
composition seemingly more convenient is possible. Namely, by means of (19) and the con- 
ditions u=c)uldn=0 on C, we have 

and hence, in consideration of (10), we obtain 

\u) 2 dxdy. 



J[u]= I (At 
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A 

Tims, taking §=§ / =L 2 (2?;l) we can decompose H into the desired form by using T c 'and 
T* c such that 

c)u 
T c u=Au, (u=^=0 onC), 

T* c v=Av. 

4.2. Before his derivation of theorem 1, one of the authors gave also the [following 
theorem 2 [10], which stands in a close connection with the former. 

Theorem 2. Let H be a self-adjoint operator of a Hilbert Space ,£), and let (a,jS) be an open 
interval containing at most a nondegenerate eigenvalue of H but no other point of its spectrum. 
Let w be any vector of the domain T) H of H with ||wj| = l, and set rj=(Hw,w), c 2 =||(fl"— y)w\\ 2 = 
\\Hw\\ 2 — rf. If e 2 <(?7— ol) (0— ??), then there is certainly an eigenvalue X' of H in (a,£) that satisfies 
the inequalities 

V- / /-^y^rj+~ (56) 

This theorem can be interpreted as follows. If we can construct a " trial function" w in 
accordance with the conditions stated above, we would obtain an estimate of the eigenvalue 
y in question by means of (56). Comparing theorem 1 with theorem 2, we note that, if a 
good choice of w is easy for the operator H under consideration, theorem 2 may be preferable 
inasmuch as it does not require the decomposition H=T*T, and also it involves fewer quan- 
tities to be calculated. However, this is not the case so far as (Pr.) is concerned, for the 
difficulties in obtaining a suitable w seem greater than those in obtaining a suitable pair of u 
and v (essentially v) because of the increased complication of the boundary conditions of 3)# . 
Furthermore, we have four components v iy v 2 , v 3 , v± at our disposal in constructing the trial 
vector function v belonging to 3) *, whereas we must make a single function w satisfy all boundary 
conditions in order that w^T) H . For instance, let us consider the case of the square plate 
dealt with in section 3 and look for a polynomial trial function w contained in QhC&I- For 
an}^ w belonging to $0?, the boundary conditions (39) of T)h reduce to 

H+^jHr ' ix==1 > - 1<y<1)> 

,2w r +45=0, (*=!, -l< y <l), 



dx 2 dif 
d 2 w 



=0' (x=y=l). 



(57) 



dxdy 
On the other hand, a general form of a polynomial of degree 2N—2, belonging to tyfl is given by 

AT m-\-n = r 

w=J2 S A mn {x 2m - 1 y 2n - 1 +x 2n -y m - 1 ). (58) 

Substituting (58) into (57), putting x=l and equating the coefficients of all powers of y to zero, 
or putting x=y=l, we derive the conditions for w given by (58) to satisfy (57). This gives 
rise to a system of 2N— 3 equations in N 2 /4 unknowns, or a system of 2N— 4 equations in 
(N 2 — 1)/4 unknowns, according as N is even or odd. The smallest value of N permitting 
nontrivial solutions is seen to be 7. Therefore, in order to obtain a function of the type (58) 
in ShOTO, we have to solve at least a system of 10 equations in 12 unknowns, leaving two 
homogeneous parameters, i. e., essentially only 1 free parameter in a 12-degree polynomial. 

4.3. Weinstein and Aronszajn have developed an approximating method [ll], 14 with a 
wide range of application, which also enables us to calculate upper and lower bounds of eigen- 

u See [11] for a standard list of references concerning this method. 
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values of a certain kind of operators. Although the arguments used in establishing the con- 
vergence of their method are theoretically interesting, it appears that, from the viewpoint of 
practical applications, the essential parts of their method are the procedures leading from 
rough bounds to those that improve as sharply as possible. In this sense their method seems 
to have some similarity to the present one, and hence it will be worthwhile to reveal the mutual 
relations between these two methods and to compare their merits. However, we reserve 
this for some future occasion and remark here only that, with respect to lower bounds of 
eigenvalues in the case of the square plate treated in section 3, the eigenvalue problem of the 
operator H can play the role of the " auxiliary problem" in their method, of which the complete 
set of the exact eigenvalues and eigenvectors is required to be known. 
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